A+/A^ Asymmetry in Hadroproduction from Heavy-Quark Recombination 
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Asymmetries in the hadroproduction of charm particles directly probe power corrections to the 
QCD factorization theorems. In this paper, the heavy-quark recombination mechanism, a power 
correction that explains charm meson asymmetries, is extended to charm baryons. In this mecha- 
nism, a light quark participates in the hard scattering that creates a charm quark and they hadronize 
together into a charm baryon. This provides a natural and economical explanation for the Aj/A~ 
asymmetries measured in nN and pN collisions. 



The production of charm particles in fixed-target 
hadroproduction experiments exhibits large asymmetries 
that are commonly referred to as the "leading particle 
effect" Q, 1^ 1^ Q . Charm hadrons that have a valence 
parton in common with the beam hadron are produced 
in greater numbers than other charm hadrons in the for- 
ward region of the beam. Asymmetries have also been 
observed in the production of light particles, such as pi- 
ous and kaons. Asymmetries in charm particles are par- 
ticularly interesting, because one can exploit the fact that 
the charm quark mass mc is much larger than Aqcd to 
make closer contact with fundamental aspects of Quan- 
tum Chromodynamics (QCD). The large mass guaran- 
tees that even at small transverse momentum the pro- 
duction process involves short-distance effects that can 
be treated using perturbative QCD. Furthermore, the 
nonperturbative long-distance effects of QCD can be or- 
ganized as an expansion in Aqcd/"^c- 

There have been many measurements of the asymme- 
tries for charm mesons yj. Several proposed charm pro- 
duction mechanisms are able to explain these asymme- 
tries by tuning nonperturbative parameters [^|^. Recent 
experiments have also measured the asymmetry for the 
charm baryon A+ |E IE Q > defined by 



a[A+] - a[A-] 



(1) 



The WA89 01 and SELEX ^ experiments observe a large 
positive asymmetry for Ac produced in the forward di- 
rection of p and beams. These asymmetries are con- 
sistent with the leading particle effect, but much larger 
than those observed for charm mesons. For tt^ beams, 
the leading particle effect predicts no Ac asymmetry, but 
a small positive asymmetry has been observed by the 
E791 and SELEX experiments. Explaining the Ac 
asymmetries is a severe challenge for any of the ^rogosed 
mechanisms for generating charm asymmetries 
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The factorization theorems of perturbative QCD [3 
imply that the cross section for A+ in a collision between 
two hadrons h,h' is given by 



da[hh' -> 



A 



(2) 



Tj/h' 



> da[ij ^ CC + X](E) 



Here is a parton distribution, da[ij cc + X] is 
the parton cross section, and D^_^^+ is the fragmen- 
tation function for a c quark hadronizing into a Ajj". 
The ellipsis represents corrections that are suppressed 
by powers of Aqcd/wc or Aqcd/p_l- The leading or- 
der processes gg — > cc and qq — > cc produce c and c 
symmetrically. Charge conjugation invariance requires 



that D„ 



= D 



SO a [Ac] = at leading order 
Next-to-leading order perturba- 



in perturbation theory. 

five corrections 0, Il0| generate asymmetries that are an 
order of magnitude too small to explain the data. There- 
fore the observed asymmetries in charm production must 
come from the power corrections to Eq. 

Recent work has shown that D meson asymmetries in 
hadroproduction and photoproduction can be explained 
by an 0(Koc,d Ijnr) power correction called heavy-quark 
recombination [ij, E E^- In the cq recombination 
process, a light antiquark q that participates in the 
hard scattering emerges with momentum of O(Aqcd) 
in the rest frame of the charm quark c and the cq 
pair then hadronizes into a D meson. In this paper, 
we extend the heavy-quark recombination mechanism to 
charm baryons. The most important process is cq recom- 
bination, which is like cq recombination except the q is 
replaced by a light quark q and the cq diquark hadronizes 
into a charm baryon. We will show that this mechanism 
can explain the observed Ac asymmetries. 
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The cq recombination cross section for a D meson is 
da[qg D] ^ da[qg cq(n) + c] p[cq{n) D] . (3) 

n 

The factor p[cq{n) D] takes into account the non- 
perturbative hadronization of a cq with color and spin 
quantum numbers n into a final state that includes the 
D meson. Since the process is inclusive, the quantum 
numbers of the cq produced in the short-distance process 
can be different from that of the D. The color and spin 
quantum numbers can both be changed by the emission 
of soft gluons in the hadronization process. The flavor 
of the recombining q can also be different from that of 
the valence antiquark in the D, but this requires cre- 
ating a light quark-antiquark pair which is suppressed 
in the large Nc limit. Neglecting such contributions, the 
heavy-quark recombination cross section for depends 
on four independent parameters: 

pi = p[cdCS^,'^) ^ D+lp, = p[cdes['^) ^ D+], (4) 
PS = pM^^f ) ^ D+],ps = p[cd{'s['^) D+]. 

Explicit expressions for these parameters in terms of 
nonperturbative QCD matrix elements can be found in 
Ref. Q|. They scale with the heavy quark mass as 
Aqcd/tIc- Analogous parameters for £)° and D~ mesons 
are obtained by using isospin symmetry and charge con- 
jugation invariance, while parameters for D*~^ states are 
related to those in Eq. Q) by heavy-quark spin symme- 
try. One might have expected the cross section in Eq. Q 
to involve a convolution with a nonperturbative function 
that depends on the fraction of the light-cone momentum 
of the D meson carried by the light antiquark q. However, 
to lowest order in Aqcd/wc, only the leading moment of 
such a distribution is relevant. Therefore, the cq recom- 
bination cross sections are calculable using perturbative 
QCD up to an overall multiplicative factor p. 

The direct cq recombination process is not expected to 
be a significant source of charm baryons, since baryon 
production requires creating at least two light quark- 
antiquark pairs and is therefore suppressed by 1 /N^ rel- 
ative to Eq. The leading recombination mechanism 
for charm baryon production is cq recombination. A lead- 
ing order Feynman diagram for this process is shown in 
Fig. n Creation of a light quark-antiquark pair is re- 
quired for the cq to hadronize into either a charm meson 
or a charm baryon, so there is a suppression in 

either case. This factor makes cq recombination a sub- 
leading mechanism for charm mesons, but the leading 
mechanism for charm baryons. The cq recombination 
cross section for Ajj" has the form 

da[qg ^ A+] = ^ da[qg cq{n) + c] ■q[cq{n) -> A+] .(5) 

n 




FIG. 1: A Feynman diagram for the heavy-quark recombina- 
tion process g + q cq(n) -|- c. The double solid lines and 
single solid lines represent charm quarks and light quarks, re- 
spectively. The shaded blob represents the hadronization of 
the diquark cq{n). 



The factor r][cq{n) A+] takes into account the non- 
perturbative hadronization of a cq with color and spin 
quantum numbers n into a final state that includes the 
A+. Isospin symmetry implies that r][cq{n) — > A+] is 
the same for q — u,d, while it is suppressed by for 
q = s. There are four possible color and spin states of 
the cq and therefore four independent ry parameters: 

= vicuCS^,'^) -> Atlm = V[cu{'s['^) A+], (6) 
V, = virM'si,'^) - A+],r)6 - v[rM's['^) ^ A+]. 

These parameters scale as Aqcd/'tIc, so cq recombina- 
tion gives a power-suppressed contribution to the cross 
section. 

The parton cross sections for cq recombination can 
be calculated using a straightforward generalization of 
the method described in Ref. |l2| for cq recombination. 
Charge conjugation is applied to the line of spinors and 
Dirac matrices associated with the recombining light 
quark in Fig. ^ Angular momentum states can then 
be projected out using the operators of Ref. 0. The 
amplitude is projected onto the appropriate color repre- 
sentation, which is either the 3 or 6 of SU{3). A simple 
prescription for projecting onto the leading moment of 
the light-cone momentum fraction of the q can be found 
m Ref. T^l. The parton cross sections for cq recombina- 
tion are 

d^ , , , ^ 2 7r^a?TO^^, 

— [qg-.cq{n)+c] = —-^-^F{n\s,t), (7) 

where s, t, and u are the standard parton Mandelstam 
variables for c + q —>■ cq{n) + c. The functions F{n\s, t ) 
for the four color and spin channels are 
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where S = s, T = t — m^, and U ~ ii ~ rn^. 

The parton cross sections for both eg and cq recom- 
bination are strongly peaked in the forward direction of 
the incoming q or q. For example, consider the ratio of 
the parton cross sections for cq recombination to that 
for gg cc, which dominates the fragmentation term in 
the cross section. We define 6 to be the angle between 
the incoming q and outgoing cq in the parton center-of- 
momentum frame. In the backward direction 9 — n, the 
ratio is suppressed by m'^/S in both ^Si channels and by 
m^/S^ in both ^S'o channels. At 6 ^ 7r/2, the ratio is 
suppressed by m^/S in all four channels. In the forward 
direction 9 = 0, there is no kinematic suppression of this 
ratio. The forward enhancement of the cq cross section 
gives charm meson asymmetries which are largest near 
xp w 1. For Ac produced mpN collisions, the fragmenta- 
tion cross section is smaller relative to cq recombination, 
so the asymmetry is large even for xp = 0.2. 

In addition to direct recombination of cq into A+, we 
need to include two additional effects: cq recombination 
into a heavier charm baryon that subsequently decays 
into and the fragmentation into A+ of a c that is 
produced in a cq or cq recombination process. The cross 
sections for A+ from the latter process, which we will call 
"opposite-side recombination" , are 

d^iig ^ + a:] = ^ da[qg cq{n) + c] 

n 

xY,p[cq{n)^D]®D^_^j^-,, (12) 
U 

da[qg ^ A+ + X] = ^ d(T[qg cq{n) + c] 

n 

xY,vmn)^B]®D^^^^. (13) 

B 

The recombination factors in Eq. (|12|l and Eq. (|13|l are 
summed over D mesons whose valence partons are cq and 
over antibaryons B whose valence partons inlcude cq. 

The feeddown from heavier charm baryons that decay 
into Ac can be taken into account through inclusive 77 
parameters defined by 

Vinc[cq{n) A+] = T^[cq{n) ^ A+] (14) 



+ ^ f^[cq{n) -> B]B[B ^ A+ + X] . 

B 

The sum over B includes all charm baryons that decay 
into A+. They include E+, S*^, and the negative-parity 
excitations A+(2593) and A+(2625) states, all of which 
have branching fractions into A+ of nearly 100%. Charm 
baryons with strangeness do not contribute to ?7inc- 

In our analysis, we choose rric = 1.5 GeV, use the 
one-loop running as with Aqcd — 200 MeV, and set 
the renormalization and factorization scales equal to 
\/ -f- m?. We use the parton distributions GRV 98 
LO for the proton and GRV-P LO 16] for the pion. 
For the fragmentation function for c A^ , we use 

^c^a+W = /a^'^(1-^) (15) 

where /^+ = 0.076 is the inclusive fragmentation prob- 
ability 17]. We also use delta-function fragmentation 
functions for the charm mesons, since this reproduces 
the shapes of their momentum distributions more accu- 
rately than Peterson fragmentation functions To'] . In the 
opposite-side cq recombination cross section, Eq. H12|) , we 
include the D and D* multiplets, but neglect the excited 
charm mesons. The best 1-parameter fit to all the D me- 
son asymmetries measured by E791 gives pi = 0.15 with 
Pi = = Ps = 0. The fit can be improved by using addi- 
tional parameters, but not dramatically. This value of pi 
is larger than the value pi = 0.06 obtained in Ref. ^3] us- 
ing Peterson fragmentation functions. Note that the sum 
of recombination parameters appearing in the opposite- 
side cq recombination cross section, Eq. \V6\ . differs from 
the inclusive parameter in Eq. (|14|l . The two are related 
if the the sum in Eq. H13|l is dominated by the lowest mass 
= and l"*^ SU(3) multiplets. Then, using charge 
conjugation and a simple quark counting argument, we 
estimate '^■B'ri[cq{n) — > S] « 3/2 r/inc[cu(ri) A+] for 
q = u,d, s. 

Because of the large uncertainties associated with the 
parton distributions of the £~, we focus on the Ac asym- 
metry data from 7r~ and p beams measured by the E791 
and SELEX experiments. The best 1-parameter fit to all 
this data yields rys.inc = 0.058, with the other three inclu- 
sive 77 parameters set to 0. The asymmetry variable a [Ac] 
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FIG. 2: The asymmetry a[Ac] for a 500 GeV n~ beam as a 
function of measured by the E791 experiment. The sohd 
curve is the best single-parameter fit with ?73,inc = 0.058, 
while the dotted curve is in the absence of cq recombination. 
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FIG. 3: The asymmetry Qf[Ac] for a 540 GeV p beam as 
a function of measured by the SELEX experiment. Fit 
parameters for the solid and dotted curves are the same as 
Fig. |5| The horizontal line at a = 1 is the physical upper 
bound. 

is shown as a function of xf for the pion beam in Fig. [3 
and for the proton beam in Fig. The 1-parameter fit 
agrees well with both the pion beam and proton beam 
data. The same fits also yield good agreement with ob- 
served p_L dependence of the asymmetries. The fits can 
be improved by using additional parameters. In Figs. El 
and 13 we also show the predictions for a [Ac] if all the r\ 
parameters are set to 0. Note that opposite-side recom- 
bination into D mesons generates a positive asymmetry 
even if all the r\ parameters vanish. It gives a reasonable 
fit to the pion beam data, but it severely underpredicts 
the asymmetry for the proton beam. Therefore, the large 
asymmetry from the proton beam is convincing evidence 
for the cq recombination mechanism. 

We have shown that heavy quark recombination pro- 
vides a natural and economical explanation of the pro- 



duction asymmetries for charm baryons as well as charm 
mesons. Further work includes a more systematic analy- 
sis of all the charm asymmetry data from hadroproduc- 
tion experiments and the prediction of charm and bottom 
asymmetries in present and future experiments. Previous 
analyses of D meson asymmetries 12, 13] do not include 
contributions from opposite-side cq recombination into 
charm baryons. This is particularly important for Ds 
mesons since any asymmetry is generated by opposite 
side recombination. 
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